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We demonstrate the existence of gravitational critical phenomena in higher dimensional electrovac 
bubble spacetimes. To this end, we study linear fluctuations about families of static, homogeneous 
spherically symmetric bubble spacetimes in Kaluza-Klein theories coupled to a Maxwell field. We 
prove that these solutions are linearly unstable and posses a unique unstable mode. We further 
I show that the associated growth rate depends only on the mass per unit length of the background 

, solution. Additionally, by a double analytical continuation this mode can be seen to correspond to 

' marginally stable stationary modes of perturbed black strings whose periods are integer multiples 

CN ■ of the Gre gory-Laflamme critical length. This allow us to rederive recent results about the behavior 

j^jQ, of the critical mass for large dimensions and to generalize them to charged black string cases. 
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I. INTRODUCTION 



^ , Bubble spacetimes were originally studied for, among other reasons, their connection with negative energy solutions 
' Qil^l^' The behavior of bubble spacetimes has been the subject of new scrutiny and use in recent years. Among 
'sJO [ these, their use in defining new solutions of Kaluza-Klcin black holes 0,0,ISI3j studies of possible phases of Kaluza- 
04 ■ Klein black holes , solutions in inflationary universe in five-dimensional spacetimes j^] and their dynamical behavior 
t ■ under different interactions via numerical simulations |^ ITol | . 

, In particular, the numerical simulations presented in Ref. 0,^3 reveal strong indications of the possible existence 
of critical phenomena in bubble spacetimes. That is, there exists a critical value A; = fee in the parameter specifying 
the strength of the U(l) gauge field in the initial data for the bubble. For this value, it was observed that, after some 
transient behavior, the solution approaches a static bubble configuration which we call the critical bubble solution. 
I I For values of k close to fee the solution approaches the critical bubble and stays close to it for some time and then 
either disperses away or collapses to a black string. The amount of proper time r that the solution remains close to 
I I the static one obeys the relation t/Rq ~ — F ln(|fc — fcc|), where Rq is the bubble location -which is fixed by the ADM 
• • . mass per unit length- and F is a universal parameter that does not depend on the initial data. 
. ^ ' present work we put on firm grounds such observation by studying linear perturbations off the critical 

k> , solutions in a g + 3 dimensional spacetime endowed with U(l) x SO(q + 1) symmetry and elucidating the existence 
5^ ' of a unique unstable exponential in time growing mode. Furthermore, we show that the growth rate fl/Ro associated 
I to this instability has the property that Q does not depend on the family of critical solutions and hence explains the 
universal behavior. We also study the limit of large q's and find that in this case the exponent f2 obeys a simple power 
law rule. 

Finally, our results can be used to draw conclusions on the 'dual black string system' obtained by appropriate 
double analytical continuations. This maps the bubble solution into a charged black string solution. Thus, 
exponential growing modes off the bubble solution turn out to correspond to static, harmonic 'deformations' of 
the uniform black strings. For these modes to exist, the strings length must be long enough to accommodate 
it, and we find that it must be an integer multiple of the critical length Lc — 27ri?o/^^- Lc results exactly that 
expected from the Gregory Laflamme instability for the uncharged black strings llj. For the charged case, however, 
the black string solutions we found are not the same as those presented in Ref. 0| as we here only consider 
U(l) gauge fields in the absence of a dilaton. Additionally, our results can be used to derive in a new way the 
exponential law for the critical mass presented by Sorkin [13 and to generalize it to the charged case. In particular, 
in the latter case we show how the addition of charge worsens the situation as far as stability of the string is concerned. 

This work is organized as follows. In section^ we summarize the obtained results and discuss the implications for 
black string spacetimes. The rest of the article presents the details of our calculations. In section UTTI we consider 
Kaluza-Klein theories coupled to a Maxwell field in spacetime dimensions q-\-3 > 5 with U(l) x SO(q + 1) symmetry. 
We perform a dimensional reduction of the action, obtaining an effective action in two spacetimes dimension, and write 



2 



down the equations of motions. In section lTVI we discuss the class of static solutions which represent the critical bubbles. 
The linear stability analysis about the bubble solutions follows in section [3 where we cast the perturbation equations 
into the form of a wave equation whose spatial operator is self-adjoint. We then show analytically that all bubble 
solutions are linearly unstable, and give analytic arguments for the existence of a unique non-degenerated eigenvalue 
A = — f2^/i?o which is negative. Furthermore, we show that for fixed spacetime dimension, the corresponding frequency 
is independent of the background solution. We provide analytical estimates for fl and compute it by numerical 
means. 



II. CRITICAL BUBBLES AND BLACK STRINGS. RESULTS AND IMPLICATIONS 



A. Critical bubbles: Summary of results 



In this article, we study the linear stability of the following static bubble solutions which are solution of the 
Einstein-Maxwell equations in {q + 3) spacetime dimensions: 



V 



2/q 



-dt^ + — + R^g 



(1) 
(2) 



where — e(l — e)(g + l)/{q^K), k denoting a coupling constant, and where 

U^l-iRo/R)''-\ V^l--e{Ro/Ry-\ 

Here, g denotes the standard metric on the g-sphere S"', and the coordinate R lies in the interval [Ro,oo). The 
parameters Rq > (corresponding to the position of the bubble) and e e [0, 1) are related to the ADM mass and to 
the asymptotic period T of the extra dimension via Madm — Rq^ r|S''^|/167r and T = 47ri?o(l — eY'^^^'^f {q — 1). 

We find spherically symmetric, in z homogeneous, linear modes that depart like exp(f2t/i?o) from the static solution. 
For the uncharged case with q = 2 this result can be drawn from the work though here we concern ourselves to 
more generic cases. For each value of the spacetime dimension parameter q we show that there is a unique such mode 
and that the corresponding dimensionless growth rate O does not depend on the parameter e which is determining 
the strength of the gauge field. This universal behavior is precisely what is expected on phy sical grounds for a 
critical solution|23- It explains the critical phenomena found in our numerical simulations P.llfl|| where the universal 
parameter F is related to Q via F = i.e. F is the inverse of the Lyapunov exponent of the single unstable 

mode. Tableland figure ^ show the obtained numerical results for as a function of q. Furthermore, our analysis 



reveals that that for q large enough, f2 is approximately given by fl 
PE^100\n-q^/'^\/q^/'^. 
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Figure |2 displays the percentile error 
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TABLE L Numerical values for Q for different values of q. 

Summarizing, bubble spacetimes admit the critical solutions I|1I2|I for spacetime dimensions q + 3 > 5 with corre- 
sponding critical exponent which is independent of e. Furthermore, for q large enough, fl ~ q^^^. Indeed, even at 
q = 7 this approximation differs only by ~ 15% from the actual critical value. 



B. Implications for black strings 

Connections between possible black string instabilities and modes in a 'dual' space obtained by Wick rotations have 
been noted, and employed, in the past (see Refs. 0,0,^3)- ^^^re illustrate how one can make the connection 
between bubble and black string spacetimes via double- Wick rotations to draw conclusions for the latter case. 

The bubble configurations (|1I2|) are related to a family of charged black string solutions via the double analytical 
continuation 1 1~> iz, z followed by the transformation e i—s- — e: 

ds' = -^dt^ + y^/" + dz^ + R^g^ , (3) 
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FIG. 1: The plot illustrates the obtained dependence of f2 on g (dark circles, with a dotted joining line for guidance). The 
inset shows versus q for low values of q. 
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FIG. 2: Percentile error PE of the actual value of Q. compared to the simple law y/q 



(4) 



where now U = 1- {Ro/R)'^-\ V ^ 1 + eiRo/R)'^-^, and = e(l + e){q + l)/{q^n). The parameters i?o (which 
now corresponds to the location of the horizon) and e satisfy i?o > 0, e > and are related to the mass and electric 
charge via 



M 



Qe 



IGtt 

{q-l)QRl-^L\S'i\ 



2e\q- 



1 



(5) 



where L denotes the asymptotic period of the extra dimension. For q ^ 2 these solutions agree with the solutions 
described by Eqs. (2.13) and (2.15) of Ref. 

Under the double analytical continuation, the unstable mode of the critical bubbles which is proportional to e^*/^° 
is mapped to a stationary mode of the corresponding black strings which is proportional to e*^^/^". This mode 
represents a static deformation of the black string, with harmonic dependency in z. That the existence of static 
modes are central in the understanding of instabilities has been discussed in |l3l , we here exploit knowing these 
modes's wavelengths explicitly to draw conclusions of black string instabilities. Clearly, for this mode to exist, the 
asymptotic periodicity of the extra dimension must be an integer multiple of Lc — 2ttRq/Q,. In the uncharged case, it 
turns out that Lc is exactly equal to the critical length for the Gregory-Laflamme instability (Compare the 

values of Table^with the values in Table 1 of Ref. j2^].) For the charged case the critical solution does not correspond 
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to the charged black string background solutions considered in[l3| as ours only include the interaction with a U{1) 
gauge field. However, as in we expect that the presence of a static deformation is related to a threshold between 
stable and unstable solutions. 

We are able to give the dependency for the critical length in terms of the parameter e in an analytic way since 
we have shown that fl is independent on e. The result is conveniently expressed in terms of the dimensionless mass 
defined in 

Li 167rlv27ry' 
For large q we can approximate fl ~ q^^^ and find 

which provides a different way of deriving the law presented in Ref. 0, |23| , and generalizes it to the electrically 
charged case. 

Note that in the charged case, the string becomes "more unstable" in the sense that the critical length at which 
instabilities arise is shorter that the charge- free case. This is to be contrasted with the case studied in IT^ where 
the opposite occurs. In that case, however, the interaction with a magnetic field was considered while in the present 
an electric field is included instead. 



2e[q-- 

q 



(6) 



III. KALUZA-KLEIN REDUCTION AND EFFECTIVE ACTION 



We consider a g + 3 dimensional spacetime (Af, g) with U(l) x SO(q + 1) symmetry. We assume that the Killing 
field dz which is the generator of the action of the U(l) group on A/ is hypersurface orthogonal. The manifold has 
the structure M — M x S'^ x with metric 



(8) 



where (p, r and gabdx°'dx^ are, respectively, a function, a positive function, and a pseudo-Riemannian metric on the 
two-dimensional manifold M. g denotes the standard metric on the g-dimensional sphere 5"^. In this article, we 
assume that q> 2. We also consider a U(l) gauge potential of the form 



Afj,dx'^ — 7 dz, 



(9) 



with 7 a function on M. 

The dynamics is described by the q 



3 dimensional Einstein-Hilbert-Maxwell action 



S^--^ f [Ric-qKg{dA,dA)]^ 
16tt J 



d" 



+3. 



(10) 



where Ric denotes the Ricci scalar belonging to the metric © and k is a coupling constant. When integrating over z, 
one obtains a reduced action describing q + 2 dimensional gravity coupled to the dilaton field and to the U(l) gauge 
potential ^ dz. Using the SO{q + 1) symmetry one can further reduce the action^^] and, after some manipulations, 
obtains 



S, 



eff 



^ g{dm,dr) 
N 



1 



q- 



r'^g{d(j>, d4)) - K r'^e-'^'l'gidj, dj) 



-gdx^ A dx^, 



(11) 



where the functions N and m are defined by = g{dr, dr), N = 1 — 2mr^ ^, and where (a;°, a;^) are local coordinates 
on M. We stress that up to this point all the fields are defined in an invariant geometrical way: e'^'^ is the norm of 
the Killing field dz, r is defined via the area radius, e~'^/'^r, of the g-spheres, and N is the norm of dr. 

In order to write down the field equations, we choose the coordinates {x'^,x^) such that x^ = r, set t — x^ and 
parametrize the two-metric according to 



dr^ 

gabdx^dx^ = -(j'^N^dt + a drf + — , 



(12) 
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where a 



with 



-g and a = —g{dt,dr)/N. The action has the form 

Seff ^-'^ j [Cc + UAdtAdr, 



a(m — am) 

1,2 



c 



M 



~2~ 



+ 



Kor 



7V(7' - ajy 



(13) 

(14) 
(15) 



where we have defined = 1/(9 + 1)/?^ a = 2q/^q + 1 for convenience, and where a dot and a prime denote 
differentiation with respect to t and r, respectively. 

Of course, we can always choose the coordinate t such that a = 0. However, we shah set a — only after varying 
the action since otherwise an equation is lost. The field equations obtained from varying H13|) with respect to a, m, 
a, and 7 are, after setting a — 0, 



2m = r'> 



N— = 



0-2 TV 



7 



2 



iV7'' 



and 



am 



Combining Eqs. pt)|) and (|17|l we also find 



1 / M / N 2m 
-{aN)' = {q-l) — 
a r'' 



(16) 

(17) 
(18) 

(19) 
(20) 

(21) 



In the next section, we discuss static solutions to these equations, where to = ■0 — 7 — 0- 



IV. STATIC SOLUTIONS 



When 7 = = we obtain m = const, so the only solutions to the field equations consists of the one-parameter 
family of g + 2 dimensional Schwarzschild spacetimes cross S^. For static configurations in general, we notice that 
Eqs. ifT^ . posses the first integrals 



(22) 
(23) 



with two constants P, Q. In fact, the static equations can be integrated completely for bubble configurations obtaining 
the following two-parameter family of solutions. 



p 

A^dx^ = —dz, 
eV 



— TT dz'^, 

y2 



(24) 
(25) 
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where = e(l — £){q + l)/{q^K), and where 

U=l-{Ro/R)''~\ V = l-e{Ro/Ry'\ 

The coordinate R lies m the open interval [i?o,oo)- The parameters Rq > (corresponding to the position of the 
bubble) and e € [0, 1) are related to the ADM mass and to the asymptotic period T of the extra dimension, see section 
nil In order to determine the latter, we replace the coordinate R by the coordinate p defined by U = and rewrite 



u 



dR' 



Therefore, in order to avoid a conical singularity at p = 0, we need T — 47ri?,o(l ~ ^Y^'^^'^ Hi ^ !)• The integration 
constants P, Q are related to Rq and P via 



P^-{q-l)Rr^P, 



Q 



{q-l)Rl 



g-1 



In the next section we consider time-dependent bubble configurations that are infinitesimally close to the static 
solutions. That is, we consider metrics of the form 



ds' = ( —e 



-2E 



1/9 



J 



where A, B, C, E are "small" functions oit and R, which vanish as i? ^ oo and such that B{t, Rq) — {\ + l/q)E{t, Rq). 
The former conditions are necessary for asymptotic flatness while the latter ensures that the period of the extra 
dimension T does not change. In order to determine the asymptotic behavior in the amplitudes a, TV, ip and 7 we 
notice that 



u 



2E 



N = gidr,dr) = U 



P , 
7 = -^ + 6r 



dr 
dR 



(26) 
(27) 

(28) 

(29) 
(30) 



Clearly, this implies that V' ~* 0, cr, — s- 1 as i? — > 00. In order to analyze the behavior near the bubble, we 
introduce the dimensionless quantity x — {R — Rq)/ Rq. The background quantities (the ones obtained by setting 
A = B = C = E = 0) have the expansion 



^a4> _ 



(i-e)2 

^2o 



0(x2 



a 

N 



Rf,''{q-l)x + Oix'), 
2q{q-iy 
(9-1)1 



2q{q-l) '^x^ [l + 0{x)], 



0(1). 



(31) 

(32) 
(33) 
(34) 



Notice that N diverges as i? Rq, but a'^N goes to zero. Linearizing about the background, we find that 6C = 
(since we do not vary r), and Sip = 2a~^6E, 5N — —2N6B. As a consequence of this, 



aSml 



^-Rl-^6B{t,RQ) = ^-^R^Q-^SE{t,RQ). 



Here and in the following, S{...) denotes the linearization of the quantity (...) about a static background. 



(35) 
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V. LINEAR PERTURBATIONS 



Since when 7 = = all solutions are static, and in view of the wave-like character of the equations H19() . H20() 
for ^p and 7, we expect that the true dynamical degrees of freedom are contained in 7 and ip. It turns out that 
this expectation is indeed true at least at the perturbative level: Our effective action Eq. H13|l satisfies the form of 
Ref. [2l| where it was shown that for linear perturbations around a static background the "gravitational" degrees of 
freedom can be completely eliminated and a pulsation equation for the perturbed matter fields alone can be obtained. 
In order to see this explicitly, we first linearize Eqs. 1)18(1 and ((16|l . obtaining 



(^iP'SiP + Ke-"'^7'(57) , (36) 

Sm' = rm {ij'Sij' + Ke""'^7'(57') + (V/^ + ne'""''-/'^) - ^ r«7Ve-"'^7'2^V- (37) 
Using r'^6N = —2r6m and the background equation a' /a — rip'"^ + Kre~"''^j''^ we can rewrite the second equation as 

[crSm]' = r'^aN {^/j'S^/j' + ne-°''''i5i) - ^ r''<7Ne-°''l'i'^5^. (38) 
Finally, using the static version of Eqs. H19() . H20|l . we find 

[aSm]' = [rViV (tA'^V + tie-°''^i 6-1)] ' . (39) 

Equations and imply that 

a5m = rVA^ + Ke-"'''7'(57) + const. , (40) 
which, using Eqs. H22|) and (|23|l . can also be written as 

aSm = K — — aPj^ 5ip + PS-f + const. (41) 

In order to determine the constant, we evaluate this equation at R — Rq. Taking into account Eq. (|35|) we obtain 
= P6^{t, Rq) + const. Since the addition of a constant to ^7 is a gauge transformation, we set the constant to zero 
in the following, which is equivalent to requiring 5j{t, Rq) = 0. Using Eq. (|41|l and the linearization of Eq. (|21|) . we 
also find 



{aNY 



(jN 



^^'^ {^'S^ + ne-'^^ySj) . (42) 



N 



Eq. (|42|l allow to reexpress all gravitational perturbations in terms of matter perturbations in the linearized versions 
of Eqs. 119(1 and 1(20(1 . Therefore, as anticipated, the linearized "gravitational" degrees of freedom can be entirely 
expressed in terms of the linearized "matter" fields, and one obtains the following pulsation equation 

Pi) - {N'^y^aNdr [aNN^PNdr {'^-^v)] +Sv = 0, (43) 



with 



1 

e-°"f' 




N = 

and where v = {Sip, ^/kS^)"^ . The spatial part of the operator acting on v in Eq. 1(43(1 is a formally self-adjoint 
elliptic operator. In the following, we start with a qualitative analysis of this operator, introducing a Hilbert space 
on which it is defined, and giving estimates for the lower limit of its spectrum. Then, we prove the existence of a 
unique bound state, i.e. an eigenfunction belonging to a negative eigenvalue. In this sense, each bubble configuration 
possesses precisely one unstable mode. The timescale associated to this mode explains the critical exponent in the 
critical dynamical behavior that has been announced in Ref. and that will be discussed in detail in Ref. [lol |. 
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A. Qualitative analysis of the pulsation operator 



Here we give a first analysis of the pulsation operator which is defined by the spatial part of the operator acting 
on V in Eq. H43|) . In order to do so, we first rewrite the pulsation equations in terms of the dimensionless quantities 
R = R/ Rf), t — t/Ro, q and e. Using the relations in section Hvl we find 



where the pulsation operator is defined by 



where 



dR ^ ' 



P Is. 



(44) 



(45) 



P = i?-? 



N 



and where 



1 
— 



2{q-lfe{l-e) 



Ri 



y2 





2(q2-l) 



9-1 



1 - 2(1 - e)y-i 



V = 1 



(l-2(l-e)y-\2V<r^) 



In order to discuss the pulsation equation, we define an appropriate Hilbert space in which the operator A (defined 
on a suitable domain) is self-adjoint. The solutions to Eq. can then be constructed using spectral theory, and 
the long-time behavior of the solutions is determined by the spectrum of A. We define the Hilbert space 

n ^ L^{{l,oo),R^;PdR) 

of all functions v : (l,oo) — > which are square integrable with respect to the measure FdR. The corresponding 
scalar product is given by 



{u,v) 



i'^Pv dR. 



Let u, V be smooth functions on (1, oo). Using integration by parts we obtain 



(m, Av) 



■^{N-^uf ■ C/N'^PN^ (N-^v) + u'^Sv 
oR oR 



dR- 



x^C/PN^ (N 

dR ^ 



(46) 



(47) 



R=l 



Defining the domain D{A) of A to be the set of all smooth functions u : (l,oo) such that all derivatives are 

uniformly bounded on (l,oo) and such that U2{R = 1) = and u and its derivatives decay faster than 
R oo, we see that the boundary term vanishes for u, v e D{A), and so 



(u, Av) ~ {Au, v) 



(48) 



for all u,v G D{A). Therefore, A : D{A) C 7i — > 7i is a densely defined, symmetric operator. Since A is real, in 
the sense that it commutes with complex conjugation, A can be extended to a self-adjoint operator; that is, A can 
be defined on a larger domain D{A) D D{A) to yield a self-adjoint operator A : D{A) C H ^ H. In general, there 
might exist different such self-adjoint extensions. If A is bounded from below, one can single out an extension by 
requiring that it possesses the same lower bound as A. In the following, we show that A is bounded from below, and 
take as extension A this unique self-adjoint extension. We refer the reader to Ref. for the theory of self-adjoint 
extensions of densely defined symmetric operators. 
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In order to show that A is bounded from below, and in order to give an estimate for this lower bound, we estimate 
the quantity 



G 



inf 

veD(A) 

^^0 



{v,Av) 
{v,v) 



(49) 



We first show that G > ~oo which implies that A is bounded from below and so there is a unique self-adjoint extension 
A which has the same lower bound G. We will also show that G is negative which means that the background solutions 
are unstable. Furthermore, we will show that A possesses precisely one eigenfunction with negative eigenvalue. In 
this subsection we provide estimates for this eigenvalue; its existence and uniqueness and the numerical computation 
of its actual value are given in the next subsection. 

Our analysis of the pulsation operator A is greatly simplified by performing the transformation v i— > Ov where the 
matrix O is given by 



O 




-2^7(1^) 
1 - 2(1 - e)V- 



(50) 



Notice that O defines an orthogonal transformation with respect to the scalar product defined in H46|) since 0"^P0 — 
P. Furthermore, O leaves the domain D{A) of A invariant. Remarkably, the transformation O brings the pulsation 
operator into diagonal form: 



where Qt = P^O^P = and 



B2 



dR dR 



dR dR 




2(g2 _ 1) 



R<i 



(51) 

(52) 
(53) 



Therefore, it is sufficient to analyze the two decoupled eigenvalue problems BiUi = XR'^ui, B2U2 = XR'^V^U~^U2, 
(mi,M2)"^ e D{A), with A e K an eigenvalue. Since B2 is a positive operator, only the operator Bi can lead to 
instabilities. Notice that Bi does not depend on e. Therefore, the bound states of the pulsation operator (if they exist) 
are independent of e. In the following, we derive lower and upper bounds for the quantity 



. /, uiBiUi dR 
mt - — - — ^ 

U3^0 J 1 J- 



(54) 



Since the transformation O is orthogonal, B2 is positive, and since we will show that Gi is negative, it follows that 
G — Gi, so Gi gives the lower limit of the spectrum of the pulsation operator. 
Using integration by parts, we find, for u e D{A): 



uiBiUi dR ~ 
1 Ji 



1 + 



dR 



> -2{q^ - 1) / R'^uldR 



which shows that Gi > —2{q^ — 1). Next, we derive upper bounds for Gi. In order to do so, we introduce the 
quantities 



E 



Z = (g-1) / RiujdR^ / {1-Uy — uldU, 



uiBiUi dR - 

JO 



dU 



l + 2±i[/ 
9-1 



dU, 
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and minimize the functional ui ^ E/Z over suitable test functions. We first study the limit q ^ oo. In this limit, E 
becomes positive. Indeed, using integration by parts and requiring Z < oo, we can rewrite 



U 



1-U d fl + U 



1 + U dU \1~U 



-ui 



dU. 



Therefore, in the limit q oo, E > and E = ii and only if ui — const • (1 — [/)/(l + U). For 3 < g < oo, this 
motivates the following test functions: 



1 - U 



(55) 



This yields 



Z = 



du < / (i-uy^du = 



and E = -l/(2g). Therefore, 



Gi<(g-lf|<- 



(g -i)(g-3) 

2q 



< 0. 



(56) 



Asymptotically as g — > cx), we have Z ~ 1/2, and the test function ui defined in Eq. ()55f) becomes optimal in this 
limit, so we expect that for q » 1, 



Gi ~ -q. 



(57) 



The numerical results in the next subsection will confirm this behavior. 

For g = 2 or g = 3, the decay of ui defined in Eq. H55|l as U approaches 1 is not rapid enough for Z to be finite. 
For this reason, we define ui = (1 — J7)™/[1 + {q+ l)U/{q — 1)] with to = 3 for g = 2 and m = 2 for g = 3. These test 
functions yield the estimates 



1 381 - 5121og(2) 
9 15-161og(2) 



-0.74 



Gi < -2 



4-31og(3) 
2 - log(3) 



-1.56 



for q = 2, 
for q ~ 3. 



(58) 
(59) 



As we will see in the next subsection, these bounds are very close to the actual eigenvalue. 

To summarize, we have shown that the pulsation operator defines a self-adjoint operator whose spectrum is contained 
in the interval [— 2(g^ — 1), oo). We have also shown that the lower limit of the spectrum is negative, and is independent 
of e. In particular, this proves that for all < e < 1, g > 2, the static bubble solutions are unstable. In the next 
subsection we show the existence of a unique non-degenerated negative eigenvalue for each q >2. 



B. Quantitative analysis — eigenvalues of A 



We consider the eigenvalue problem 



- {UR'- 



w 



2(g2 _ 1) 



1 



XR'^w, 



(60) 



where here and in the following, a prime denotes differentiation with respect to R, and w is a smooth function on 
(l,oo) such that all its derivatives are bounded near ^ = 1 and decay faster than ^""J/^-i/^ near R = oo. We 
show that there exists precisely one negative eigenvalue A to this problem and compute it by numerical means. The 
following analysis is standard and follows the theory of singular Sturm-Liouville problems |23l |. 

We first show that there exists local solutions with the required behavior near ^ = 1 and R = oo. At the location 
of the bubble {R = 1), we have a regular singular point while at = oo we have an essential singularity. In order to 
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analyze the behavior near i? = 1, it is convenient to replace the coordinate Rhy U = 1 — R ^\ We then obtain 
the system 




Therefore, to leading order in U, Eq. (|61|l reads 

9$ 



d{\og U) 



with solution $0 = (1 + log [/)eo, where cq is a constant vector. Since we require $ to be bounded near [/ = we 
have to choose Cq = (1, 0) in order to avoid the logarithmic term. According to the general theory of regular singular 
points [23|j this gives rise to a one-parameter family of local solutions near ?7 = of the form $ = '^IkXd^kU^ ■ 
Plugging this form into Eq. H61|l yields the recursion relation 

k 

efc = (fcl- Ao)-i^A;efc_i , /c= 1,2,3,... , (62) 

1=1 

where I denotes the identity matrix. The series converges for t/ > small enough |23| . 

In order to analyze the behavior near R = 00, we first introduce the coordinate p defined by 

which runs from zero to infinity. Notice that p/R 1 as R 00. Therefore, there is a smooth function h{p) such 
that h{p) — > 1 as p ^ 00 and such that p = h{p)R. In terms of this new coordinate, we can rewrite the eigenvalue 
problem as 

—X + {X + B{p))X^O, B{p) = --— + ^ '—^ (64) 



R'^1 



1 + 



where X = Tw,T^ R'^/^U^'^. Since B{p) — 0{p ^) it follows from a theorem by Dunkel ^23\ that there exists local 
solutions near p = 00 of the form 

X(p) = e-^qi + c(p)], (65) 

where VL — V— A and where c{p) is a twice continuously differentiable function which satisfies c{p) and dpc{p) 
as /9 — > 00. The introduction of the "tortoise" coordinate p is important for the case q = 2, where 



P 



^Vr-i + log (^Vr + Vr-i 



and so e~^P — j^-^/'^^-^R \p 4- c{R)\ , where h = 2^^e^/^ and c{R) is twice continuously differentiable and satisfies 
c{R) 0, dpjc{R) — > as /J ^ 00. For q > 3 one can show that h = \ + 0{R^^), and in those cases, one can simply 
replace e'^^ by e~^^ in the expression H65|l . 
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Therefore, there exists local solutions wl^wh, near left and right boundary points respectively, of the form 

WLiX; R) = 1- (^2i±j + (^-^) U + 0(C/2), (66) 
wr{X;R) - U-^/^R~'^/^e~^P[l + c{X-p)], (67) 
where c(A; p) is a twice continuously differentiable function which satisfies c(A; p) — s- and dpc{X] p) — > as p — s- cx). 



By varying A in the range of negative real values, we can look for a global solution by the "shooting and matching" 
procedure. 

In what follows, we will prove that there exists such a global solution and that it is unique up to a multiplicative 
factor. The idea is to use the fact that the local solution at ^ = 1 for A = 0, which is given explicitly by 

wl{0; R) = wo{R) ^ lZ\l[ll^^l , (68) 

has exactly one node in the interval < ^ < oo and use the nodal theorem which states that the number of bound 
states is equal to the number of nodes of the local solution belonging to zero eigenvalue. For a precise statement of 
this theorem, which also applies to coupled systems of Schrodinger equations, we refer to p^ . Below, we provide an 
elementary proof for the existence of a unique negative eigenvalue. This proof consists of three parts: (i) To show 
that the eigenvalues are non-degenerate assuming they exist; (ii) that an eigenvalue exists, by demonstrating that 
the 'shooting and matching' procedure has a solution and (iii) that the eigenvalue is unique. 

In order to proof these assertions we rely on the following observations. Let wi and W2 be two local solutions to the 
eigenvalue problem (|6U|I . belonging to the eigenvalues Ai and A2, respectively. The Wronski determinant is defined by 

W[wi,W2] = UR'^ {w[w2-w'2Wi) , (69) 

and by virtue of Eq. (|60|l . it satisfies 

W'[wi,W2] - (A2 - Xi)R''wiW2 . (70) 
The Wronski determinant will play a central role in proving the three statements above. 

(i) The eigenvalues are non-degenerate. Suppose there are two eigenfunctions, wi, W2 belonging to the same 
eigenvalue Ai — A2. Then, Eq. H70I) implies that W^[wi,W2] = const. But this constant is zero since wi, W2 and 
their derivatives are bounded near R — 1, and because of the factor U in the definition of the Wronski determinant. 
Therefore, Vl^[?«i,W2] = 0, which means that wi and W2 are linearly dependent. 

(ii) There exists a negative eigenvalue. Let ^1 > 1 be large enough such that for each —2{q^ — 1) — 1<A<0 and 
all R > Ri, wb{X;R) > 0, w'j^{X;R) < and wo{R) < 0. That we can always choose such a ^1 follows from the 
definition of wq and the behavior of wr near the right boundary. Additionally, we can extend wl{X; R) to R = Ri. 
Now define W{X) = M^[wl(A), Wfl(A)]|j^^^^. We can match wl(A) and wr{X) if and only if W{X) — 0. So our task 
it to show that the function A 1-^ W{X), which is continuous, possesses a zero. Now choose A = Ai = —2{q'^ — 1) — 1. 
Since 

WLiXi;R) = l+ . ^ +0{U^), (71) 
W - 1)^ 

we have wl{Xi; 1) = 1, w^(Ai; 1) > 0. Furthermore, it follows from Eq. that for A < — 2(g^ — 1), w{R) cannot have 
a local maxima as long as w > 0. Therefore, wl{Xi;R) is monotonously increasing and we must have wl{Xi;Ri) > 0, 
w^(Ai;^i) > 0. This implies that Vl^(Ai) > 0. On the other hand, for A = 0, we have wl{0;Ri) < 0, w'j^{0\ Ri) < 
which implies W{Xi) < 0. Therefore, there exists A with Ai < A < for which W{X) ~ and for which the lo- 
cal solutions wl(A) and wr{X) can be matched to yield an eigenfunction wx{R) belonging to the negative eigenvalue A. 

(iii) The negative eigenvalue is unique. In order to prove this, we first show that w\{R) possesses no nodes in 
(1, 00). We normalize w\ such that W\{1) = 1. Suppose that wx has a node at .R2 > Ij and that w\ is strictly positive 
on (1, ^2)- Then, u'^(^2) < 0. Integrating Eq. H7()(l from _R = 1 to = ^2 we have, choosing Ai = A and A2 = 0, 

[C/i?«w>o]^^^^ = W[wx,wo]{R2) - W[wx,wo]{l) ^ -X J R^w^wodR. (72) 
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It follows that 'Wo{R2) < for otherwise, the left-hand side would be negative or zero and the right-hand side positive. 
Therefore, lies to the right of the single node of wq. Next, integrate Eq. (|70|l from R = R2 to some R > R2 and 
obtain 

/•« _ 

URiw'^{R)wo[R) = "A / R''wxWodR+ [C/i?«<wo] ^^^^ + UR'iwx{R)w'o{R). (73) 

Suppose first that w\ <Q on (i?2, R)- Then, the right-hand side is positive, so we must have w'^{R) < 0. On the other 
hand, if w\ > on (i?2, -R) the right-hand side is negative, and we must have w'^lR) > 0. Both options contradict the 
fact that 'W\{R) — > as ^ ^ 00. Therefore, wx is positive on (l,oo) and possesses no nodes. Finally, suppose now 
there exists another eigenfunction belonging to a negative eigenvalue A > A. Integrating Eq. (|70|l over (l,oo) and 
taking into account the boundary conditions, it follows choosing Ai = A and A2 = A in Eq. H70(l . that 

/oo 
WwxW'^dR, (74) 

which implies that Wj^ must have at least one node. But then we can repeat the argument in the first part of (iii) and 
reach a contradiction by replacing wx by w^. 

Therefore, there exists a unique non-degenerated negative eigenvalue. Although we have not proven that this 
eigenvalue corresponds to the ground state, we conjecture that the spectrum of the pulsation operator consists of the 
single negative eigenvalue and a continuous spectrum [0, cxd). 

We have computed the eigenvalue numerically following the "shooting and matching" procedure described in |27l |. 
In order to do so, we rewrote equation Ht)4|l in terms of the variable Y = e^^X and implemented it numerically. 
To proceed with the numerical integration we introduce the compactified coordinate U = 1 — R~'^'^~^^ and define 
left /right boundaries at ?7 = + A and J7 = 1 — A respectively. At these boundaries, initial values are given in 
correspondence with the local solutions H66I) . H()7|) . respectively, where the left local solution is truncated after the 
second order term in U and where c{p) = is chosen in the right solution. The equations are then integrated using 
the LSODE package combined with a nonlinear Newton iteration to match the solutions at an intermediate point 
by modifying the parameter A and the scaling of the right solution. The tolerance specified in the LSODE solver, 
which is used in the variable time-step ODE integrator, was set to TOL = 10~^^. Additionally, we require that at 
the intermediate point, Y and its derivative match with an error smaller than 10~^. 

As an illustration, figure |3 presents the obtained solutions Y = e^^'^X for the cases q = 2,7, 20. The matching point 
was chosen sd U = 0.6 but we checked that the obtained results are independent of this choice. 

We then use the developed code to obtain the dependence of fl on the dimensionality q in the range q = 2, 3, ...250. 
To ensure convergence of the obtained value we run the code with A = 0.01/2" (n an integer number) and increase n 
until the difference of two successive values of 57 lies below 0.01%. The results for are presented in figure while 
figured displays the percentile error PE — 100 \ y/q— i}\/^/q. As expected, the discrepancy diminishes as q increases. 



VI. FINAL WORDS 



This work has put on firm grounds the existence of critical phenomena in bubble spacetimes confirming the indica- 
tions obtained via numerical simulations. We find spherically symmetric, in z homogeneous, linear modes that depart 
like exp{nt / Rq) from the critical solution which is static's^]. We have identified the critical solution, the growth rate 
of the unstable mode and established the universality of the growth rate. The universality follows from the fact that 
it is possible to decouple the pulsation equations into two Sturm-Liouville problems, where one problem does not give 
rise to instabilities while the other one is manifestly independent of the background solution when expressed in terms 
of appropriate dimensionless quantities. Additionally, the study of the unstable Sturm-Liouville problem allows for a 
derivation of the simple relation Q ~ q^^^ in the limit of large g's. This relation agrees very well with the numerical 
results. 

Furthermore, we have made use of a double analytical continuation to map the solution to a stationary mode of the 
corresponding black strings which is proportional to e*^^/^" . This mode represents a static deformation of the black 
string, with harmonic dependency in z. This mode is allowed if the asymptotic periodicity of the extra dimension is 
an integer multiple of Lc = 2ttRq/^. For the uncharged case, the obtained values of give the the critical length 
for the Grcgory-Laflamme instability In the charged case, the black string solutions given by H1I2I) do not 
correspond to those considered in 12] since we restricted ourselves to a U(l) gauge field in the absence of a dilaton. 
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FIG. 3: The eigenfunction Y = e^^X as a function of U for three different values of q. 



We give the dependency for the critical length in terms of the parameter e in an analytic way. 
expressed in terms of the dimensionless mass introduced in[r^ 



IGtt 



9-1 



2e[q-- 

q 



This result can be 



(75) 



For large q we can approximate f2 ~ q^^^ and find 



l + 2e), 



(76) 



16 \2ttJ 

which provides a different way of deriving the law presented in Ref. 0,1131, and generalizes it to the charged case. 

Note there is a remarkable difference between the instability of the charged strings considered here and those in 
[T^ll^ . Whereas the addition of charge makes the critical length, at which the instability arises, shorter in our case 
-making the string "more unstable"- while longer (and the extremal case marginall y st able) in 12, 19] . The reason 
for such a marked different is likely due to the interaction with a magnetic field in |12. [l9l | while an electric field in 
the present case. Thus, the dynamics of charged strings even in 5-dimensional cases appears to be quite rich. 

Finally, (single) Wick rotations are often used to study and understand phenomena in the Euclidean sector and 
draw conclusions about black branes. Here, and also in it is show that the understanding of solutions of bubble 
spacetimes can be used to shed light into black string systems and viceversa. 
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